Mon. Not. R. Astron. Soc. 000, HHTOl (20111 Printed 21 January 2013 (MN WT&t style file v2.2) 



A cosmic speed-trap: a gravity-independent test of cosmic 
acceleration using baryon acoustic oscillations 
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ABSTRACT 

We propose a new and highly model-independent test of cosmic acceleration by com- 
paring observations of the baryon acoustic oscillation (BAO) scale at low and interme- 
diate redshifts: we derive a new inequality relating BAO observables at two distinct 
rcdshifts, which must be satisfied for any reasonable homogeneous non-accelerating 
model, but is violated by models similar to ACDM, due to acceleration in the recent 
past. This test is fully independent of the theory of gravity (GR or otherwise), the 
Friedmann equations, CMB and supernova observations: the test assumes only the 
Cosmological Principle, and that the length-scale of the BAO feature is fixed in co- 
moving coordinates. Given realistic medium-term observations from BOSS, this test 
is expected to exclude all homogeneous non-accelerating models at ~ 4cr significance, 
and can reach > 7a with next-generation surveys. 

Key words: cosmology - dark energy 



1 INTRODUCTION 

In the last 10-15 years, the ACDM model has been es- 
tablished as the standard model of large-scale cosmology; 
the model is an excellent match to many observations in- 
cluding the anisotrop ies in the CMB measured by WMAP 
|Komatsu et aj|201ll ) and other experimen ts, the large-scale 
clustering of galaxies |Percival et alll2~010ft, th e Hubble dia- 
gram for high-z supernovae ( Guv et all l20ld ; IConlev et all 
l201ll ). and the abundance and b aryon fraction of rich clus- 
ters of galaxies (|Allen et alll201ll ). 

Despite these great observational successes, the model 
appears unnatural since 96% of the universe's mass- 
energy is not observed, but is only inferred from fitting 
the observations. Also, the dark sector contains at least 
two apparently unrelated components, dark matter and 
dark energy; recent review s of d ark energy a re giv en by 
iFrieman. Turner fc Hutererl (120081 ) and lLinderl (|2008ft . 

The most direct evidence for cosmic accelerat ion comes 
from the Hubble diagra m of Type-la supernovae |Guv et all 
l20ld ; IConlev et al201lft . which shows that SNe at 0.3 < z < 
0.9 are fainter, relative to local SNe, than can be accommo- 
dated in any Friedmann-Robertson- Walker model without 
dark ene rgy. A model-independen t approach has also been 
given bv lShapiro fc Turner! (|2006t ). who show that the SNe 
results require accelerated expansion at z < 0.4 at around 
the 5a significance level without assuming the Friedmann 
equations. 
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However, there are some possible loopholes in the super- 
nova results: since they are fundamentally based on bright- 
ness measurements, the interpretation could be affected by 
either unexpected evolution of the mean SNe properties over 
cosmic time, or some process which removes photons en 
route to our telescopes, such as peculiar dust or more ex- 
otic effects such as photon-dark matter interactions. The 
simplest such effects with monotonic time-depe ndence are 
stron gly disfavoured by SN observations at z > 1 (Ric ss et all 
12007ft . but more complex time-dependent effects could still 
leave these loopholes open. 

Independent of supernovae, there is powerful sup- 
port for dark energy from observations of the anisotropies 
in the cosmic micr owave background l|Larson et a il I2010I : 
iKomatsu et alll201lft an d the large-scale clustering of galax- 
ies (|Percival et alll2010ft . but this is dependent on assuming 
general relativity and the Friedmann equations; if these both 
hold, the model parameters are tightly constrained by CMB 
and LSS data, and the expansion history a(t) must match 
ACDM models within a few percent. However, in alterna- 
tive gravity theories, we cannot make model-independent 
statements from the CMB or large-scale structure: clearly 
any successful modified-gravity model should eventually be 
consistent with these observations, but the model space of 
modified gravity is large and the calculations non-trivial; 
so in non-GR models we cannot necessarily use the CMB 
and LSS observations to make any definite statement about 
recent acceleration. 

The accelerated expansion is so startling that it is de- 
sirable to test it via multiple routes with a minimum num- 
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ber of model assumptions. A very direct test of accelera- 
tion has been proposed using the "cosmic drift", which is 
the s mall change in r edshift for fixed object(s) over time 
(e.g. iLiske et~al l200si ); the predicted change is dz/dt = 
(1 + z)Ho — H(z). However, this effect is tiny over human 
timescales, of order cm/s/year, and will probably require 
over 20 years baseline to get a convincing detection. 

Here we propose a new and robust test for cosmic ac- 
celeration based only on the cosmic "standard ruler" in 
the galaxy correlation function: in the standard model, this 
is a feature created by acoustic oscill ations in the baryon - 
photon fluid before recombination (e.g. [Peebles fc Yu|[l970j); 
this w as analysed in more detail by lEistenstein fc Hul 
l|l998l) and lMeiksin White fc Peacockl (ll999T) . then first de- 
tecte d in 2005 by lEisenstein et all (|2005h in SDSS data, 
and ICole et all (120051 ) using the 2dFGRS survey. The 
length of this ruler, hereafter r„, depends only on mat- 
ter and radiation densities and is accur ately predicted 
from CMB observations at « 153 ± 2Mpc |Komatsu et all 



. Many recent studies (e .g. Eiscnstcin, Sco & Whit' 
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Tian et al 201 lh have shown how precision measurements of 
this BAO scale from huge galaxy redshift surveys can pro- 
vide powerful constraints on the properties of dark energy, 
and test for evolution of dark energy density; more details 
are given in Section [2] 

However, in the current paper we do not assume any 
gravity theory or the actual length scale of this feature, only 
that we can observe some feature at a specific lengthscale 
imprinted on the galaxy distribution at high redshift, which 
expands with the Hubble expansion and remains a constant 
ruler in comoving coordinates. We then derive an inequal- 
ity relating observations comparing this ruler at low and 
intermediate redshift, which is satisfied in any reasonable 
non-accelerating model, but is violated by accelerating mod- 
els approximating ACDM. In more detail, we use the radial 
component of the BAO feature at Z2 ~ 0.75 to constrain the 
product H(z2)r a , and we then compare to the spherical- 
averaged BAO feature at low redshift z\ ~ 0.2, which is 
related to the average of 1/H(z) at ^ z ^ Z\. Then, as- 
suming any non-accelerating model we derive a strict upper 
limit on the ratio of these. Models approximating standard 
ACDM predict a result which violates this inequality by a 
substantial amount ~ 10 — 20%, depending on cosmologi- 
cal parameters and redshift. Future large redshift surveys 
should be able to measure this ratio to ^ 2% precision: as- 
suming our inequality is significantly violated as predicted, 
we can then exclude all homogeneous non-accelerating mod- 
els regardless of Friedmann equations, gravity theory or de- 
tails of the expansion history. 

The plan of the paper is as follows: in § [2] we review 
the basic features and observables of baryon acoustic oscil- 
lations. In § [3] we derive the new inequality relating BAO 
observables for non-accelerating models. In § [4] we discuss 
future observations and related issues, and we summarise 
our conclusions in § [5] 



2 OBSERVATIONS OF THE BAO FEATURE 

The baryon acoustic oscillation (hereafter BAO) featu re 
jEistenstein fc Hulll998l : lMeiksin. White fc Peacockl! 19991 ) is 



a bump in the galaxy correlation function £(r), or equiva- 
lently a decaying series of wiggles in the power spectrum 
P(k), corresponding to a comoving length denoted by r s , 
created by acou stic waves in the early universe prior to 
decoupling. (See iBassett fc Hloze"kl l|2010l ) for a recent re- 
view). In the standard model, its length-scale is essentially 
set by the distance that a sound wave can propagate prior 
to the "drag epoch" at Zd ~ 1020, denoted r s (zd), and 
this length depends only on physical densities of matter 
f2 m /i 2 and baryons Q,bh 2 , (together with radiation density 
Q r h 2 which is pinned very precisely by the CMB temper- 
ature). In the standard model the relative heights of the 
acoustic p eaks in CMB anisot ropics constrain Q. m h 2 and 
Qbh 2 well |Komatsu et a3 l201ll ). which leads to a prediction 
r s ~ 153 Mpc comoving with approximately 1.5 percent pre- 
cision. This predicted length does not rely on the assumption 
of a flat universe, since the relative CMB peak heights con- 
strain the various densities reasonably well without assum- 
ing flatness. However, the CMB-predicted length rs does 
depend on assuming standard GR, and several assumptions 
about the mass-energy budget including standard neutrino 
content, negligible early dark energy, no late-decaying dark 
matter, negligible admixture of isocurvature perturbations, 
etc. However, in the rest of this paper we leave r s as an 
arbitrary comoving scale, which cancels later. 

The BAO feature provides a standard ruler which can 
be observed at low to moderate redshift using very large 
galaxy redshift surveys; in the small angle approximation 
and assuming we observe a redshift shell which is thin com- 
pared with its mean redshift z, there are two primary ob- 
servables derived from a BAO survey: firstly the angle on 
the sky subtended by the BAO feature transverse to the 
line of sight, A8(z) — r s /[(l + z)Da(z)], where Da(z) is the 
conventional (proper) angular-diameter distance to redshift 
z; and secondly the difference in redshift along one BAO 
length along the line of sight is Az//(z) = r„H (z)/c (e.g. 
iBlake fc Glazebrooklliooll . ISeo fc Eisensteinll2003l ). We note 
that calculating comoving galaxy separations from observed 
positions and redshifts requires a reference cosmology, hence 
a difference between the true and reference cosmology will 
produce an error in the inferred r s ; however, any error in 
the reference model cancels to first order in the dimension- 
less ratios r„/ ' Da{z) and r s H(z)/c, so both of these ratios 
can be well constrained with minimal theory-dependence by 
measuring BAOs in a galaxy redshift survey. 

The ability to independently probe Da(z) and H(z) is a 
powerful advantage of BAOs over other low-redshift cosmo- 
logical tests. Furthermore, a redshift survey useful for BAOs 
can also measure growth of structure via redshift-space dis- 
tortions and thus test for consistency with GR, though we 
do not consider this here. 

However, in practice, current galaxy redshift sur- 
veys are not quite large enough to robustly measure the 
BAO feature separately in angular and radial directions 
(thou gh there are tentative detections, e.g. iGaztanaga et all 
120091 ). The current measurements primarily constrain a 
spherically-averaged s cale, called Dv, which is defined by 
lEisenstein et all (|2005l ) as 



Dv(z) = 



{l + z) 2 D 2 A {z '"■ 



1 1/3 



H(z) 



(1) 
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this is essentially a geometric mean of two transverse di- 
rections and one radial direction. Observations using the 
2dFGRS and SDSS-II redshift surveys have measured the 
dimensionless ratio d(z) = r B /Dy{ z) at low redshifts 
(IPercival et a 3 l20ld : iKazin et all |2010| ). which we discuss 
later. We note that as z — > 0, Dy(z) — > cz/ Hq; however, this 
approximation is not very useful in practice, since we can- 
not measure the BAO feature at very low redshift z < 0.02 
where corrections of order z 2 axe unimportant. We give a 
better approximation below in § 14.21 

In practice, the BAO feature is not a sharp spike but 
a hump in £(r) of width approximately 15% of r s , so there 
are several subtle effects in actually extracting the scale r s 
from a redshift survey: we discuss these in more detail in 
§ 14.11 However, for the purposes of this paper we only need 
to assume that r s is a constant comoving length to ~ 1% 
at redshift ^ 0.8, so these precision details are relatively 
unimportant for the rest of this paper. 



3 THE COSMIC SPEED TRAP 

Here we derive a new inequality which we denote the "cos- 
mic speed-trap", which must be satisfied by any reason- 
able non-accelerating model, but is violated by ACDM and 
other accelerating models. We start off by assuming an arbi- 
trary non-accelerating model, and deriving a lower limit for 
Dy{zi) in terms of the value of H(z 2 ) at a higher redshift z 2 . 
Then, we form a ratio of BAO observables which eliminates 
H(z2) and r B , and we obtain the speed-trap inequality (fl"5|) 
which forms our main new result. 

3.1 An inequality for Dy in non-accelerating 
models 

Here we derive an inequality for Dy(z) which is satisfied in 
any non-accelerating model, but may be violated by accel- 
eration. 

First we define as usual a to be the cosmic expansion 
factor relative to the present day with ao = 1, redshift z by 
1 + z = a -1 , and the Hubble parameter H(a) — a/a where 
dot represents time derivative. Then we have the expansion 
rate 

a = aH(a) = . (2) 

1 + z 

if the expansion of the universe was non-accelerating, then 
a is non-positive and the function above must be non- 
increasing with time or a, therefore non-decreasing with 
increasing z. Therefore, if we consider any two redshifts 
z 1 < z 2 , in any non-accelerating universe, 

H(zi) H(z 2 ) 

Assuming only the cosmological principle, any observed vio- 
lation of this inequality is a direct proof that the expansion 
has accelerated, on average, between the earlier epoch z 2 and 
the later epoch zi, without reference to any specific theory 
of gravity or geometry. 

A concordance ACDM model does violate this inequal- 
ity due to the recent positive acceleration: a minimum value 
of H(z)/(1 + z) occurred at z acc = y/2 Q,\/Q, m — 1; for 
the concordance value Q m ~ 0.27, this gives z acc « 0.75, 




0.2 0.4 0.6 0.8 1 



Figure 1. The expansion rate a = H(z)/(l+z) from the standard 
Friedmann equation is shown as a function of redshift, for four 
models: three fiat ACDM with respectively Q m = 0.24 (dashed), 
0.27 (solid) and 0.31 (dot-dash); also one wCDM model with 
f^m = 0.32, w = —0.85 (dotted). The upper panel shows expan- 
sion rate relative to present-day Ho, while the lower panel shows 
absolute expansion rate, where for each model Ho is adjusted to 
keep the CMB acoustic angle constant. Parameters are given in 
Table [T] 

and H(z iscc )/(1 + z aC c) w 0.85 Ho- The expansion rate 
H(z)/(1 + z) is shown in Figure [T] for a few representative 
models: it is notable that the value of H(z)/(1 + z) remains 
within a few percent of its minimum between 0.5 ^ z ^ 1.2, 
and it rises rather sharply at low redshift; for the concor- 
dance model it only crosses the half-way value between the 
minimum and the present-day Ho at the modest redshift of 
z w 0.17, and three-quarters of the speedup has occurred 
since z w 0.31. Thus the actual speedup of the expansion 
rate is quite concentrated at rather low redshift; this be- 
comes relevant later. 

Next, we suppose we have a measurement of H(z2) at 
an earlier epoch z 2 ; for a non-accelerating model we now 
derive a lower limit on Dy(zi) at a later epoch z\ where 
Zi < z 2 . 

The comoving radial distance to redshift zi is 
D^)=c£ 1 ^ )d z • (4) 
If the universe is non-accelerating and zi < z 2 , we can rear- 
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Model 




H 


w 


to 






(kms -1 Mpc" 1 ) 




(Gyr) 


C 


0.27 


70.0 


-1 


13.86 


L 


0.24 


72.5 


-1 


13.82 


H 


0.31 


67.1 


-1 


13.91 


W 


0.32 


64.6 


-0.85 


13.98 



Table 1. Cosmological parameters for the four example models 
discussed in the text; model C is the baseline concordance model, 
while the others are selected to roughly span the current 2a al- 
lowed range in Q m and w. All are flat, and have Hq adjusted to 
give very similar values of ij\ consistent with WMAP, therefore 
have similar values of to. 



range inequality © into 1/H(z) > (1 + z 2 )/[H(z 2 )(l + z) 
inserting this we have 



H(z 2 ) 



(•») 



The proper angular-diameter distance Da{z) is defined 



by 



(1 + z)D A (z) S \Ro\ S k J = D R (z)^l (6) 

where \Rc\ is the curvature radius of the universe in co- 
moving Mpc, x = Dr(z) /\Rc\, and the function Sk(x) = 
sinhx, x, sin a; for the cases k = — 1, 0, +1 where k is the 
sign of the curvature. 

Note that in the above we have left Rc as a con- 
stant but arbitrary curvature radius, thus we have not 
assumed the Friedm ann equation which gives Rc = 
(c/Ho)y // k/({l to t — 1); we have only assumed that the uni- 
verse has a metric with some well-defined curvature ra- 
dius Rc, whic h follows from t he assumption of homogeneity 
and isotropy (|Peacockl I1999T ). Also, we have not assumed 
any functional form for H(z), only that it obeys the non- 
acceleration condition ((3| at all z ^ z 2 ; what happened ear- 
lier at z z 2 is immaterial. 

For the other term in Dv, we use a similar inequality 
for 1/H(zi) as above, which is 



CZ\ 



C2l(l + Z 2 ) 



(<) 



H(zi) ' H{z 2 )(l + Zl ) ' 

substituting both of the above into Eq. [TJ we obtain the 
inequality 

1/3 /„ , x \ 2/3 



Dv(zi) > 



c(l + z 2 ) 
H(z 2 ) 



Zi (ln(l + 2i)) 2 



1+Zl 



Xl 



.(8) 



where Xi = Dr(zi)/\Rc\ as above. 

This inequality is strict for any non-accelerating and 
homogeneous universe with a Robertson- Walker metric, in- 
dependent of details of the expansion history or the gravity 
model. This is not so useful on its own, but we will see in 
the next section how to combine observables to cancel the 
z 2 dependence. 

We note that the factor (Sk(xi) /xi) 2 ^ 3 = 1 exactly for 
flat models, and is ^ 1 for open models (so open models 
always strengthen the inequality); the factor is ^ 1 for closed 
models which weakens our inequality, but only by a small 
amount if we consider sufficiently low redshift Zi, since the 



effect of curvature on distances only enters to third order in 
z; at small x and k — +1 we have 



Sk(x) 



2/3 



1 - 



(9) 



therefore we need an upper limit on x for closed models. We 
get a firm limit as follows, using an upper bound on Dr and 
a lower bound on Rc for closed models. 

To limit Dr, we can use the non-acceleration inequality 
© between z = and an upper redshift z 1 to get 1/H(z) ^ 
l/[Ho(l + z)], which now leads to an tipper bound on Dr(zi) 
in terms of Ho, Dr(zi) ^ (c/f/o)hi(l + Zi) ^ czi/Ho, for 
any non-accelerating model. This gives x\ ^ czi/ HoRc ■ 

We may also obtain a lower bound on Rc as follows: in a 
closed model, it is clear from Eq.[U]and sin a; ^ 1 that Da(z) 
cannot exceed Rc /(1 + z) regardless of the expansion history 
H(z). If we take for example Rc = 0.6 c/ Ho and Ho = 
70 kms" 1 Mpc" 1 , this leads to D A (z = 3) < 642 Mpc, only 
0.4x the concordance value of 1638 Mpc. However, observed 
angular sizes of z ~ 3 galaxies already convert to rather 
small physical sizes based on the concordance model, and 
making them smaller by another factor < 0.4 appears to 
be seriously discrepant. We therefore exclude closed models 
with R c < 0.6 c/ H . 

A stronger lower bound may be obtained with other 
methods: e.g. t he luminosity dis tance Dl(z = 1.5) mea- 
sured from SNe l|Riess et "al 120071 ) agrees well with the con- 
cordance model, and if we adopt a lower bound 0.8 x the 
concordance value, we obtain Rc ^ 0.84c/-ffo- However, to 
remain fully independent of SNe data we do not use this 
below. A stronger limit should also be possible in future 
using angular BAO measurements at z ~ 3, e.g. from the 
HETDEX or BOSS projects. 

However, for the following we take Rc ^ 0.6 c/Ho as a 
conservative gravity-independent lower limit for closed mod- 
els. This leads to a firm upper limit x ^ ln(l + z)/0.6 for 
closed non-accelerating models, which we use below. 

3.2 The observable speed-trap 

The above inequality ((8} relates the volume-distance Dv{zi) 
at low redshift to the Hubble constant H(z 2 ) at a higher 
redshift. Neither of these quantities are directly observable 
at present, but it is possible to measure both of them relative 
to the BAO length-scale r s ; then, dividing these two cancels 
the length scale r s and gives a ratio measurement. Applying 
the Dv inequality above gives us a limit which must be 
satisfied by any reasonable non-accelerating model, but is 
found to be violated by an expansion history close to ACDM, 
for a range of suitable choices of z\ ~ 0.2, z 2 ~ 0.75. 

The Hubble parameter H(z 2 ) may be measured using 
the radial BAO scale (along the line of sight) in a redshift 
shell near z 2 ; for a thin shell and ignoring redshift-space 
distortion effects, this gives the observable 



Az//(z 2 ) 



r s H(z 2 ) 



(10) 



In practice it is useful to divide by 1 + z 2 and define 

- - ^\ (ID 

1 + 22 C(l + Z 2 ) 

since this y is rather close to a constant over a substantial 
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range of redshift in a ACDM model (as in Figure [T]), and we 

will see that it has a convenient cancel lation below. 

Using the SDSS-II redshift survey. iPercival et all (|2010l ) 
have already measured the dimensionless ratio 



d(z) = r s /D v {z) 



(12) 



at redshift z — 0.2 and 0.35, and also a combined ratio at 
z — 0.275. (We discuss the numerical results later). 

We now form the ratio of observables zid(zi)/y(z 2 ) 
which gives, from the definitions above 



"(I + 22) 



Z\ 



(13) 



zid(zi) 
V{Z2) H(z 2 ) D v {zi) 

assuming only that r s is a fixed comoving ruler independent 
of z. 

If we now assume that the universe has never acceler- 
ated below redshift 22, we may apply the inequality ((8J for 
Dy(zi); this cancels the 22 factors, giving the inequality 

i/3 / \ 2/3 

(14) 



Zi d(2i) 

viz*) 



zt{l + z 1 ) 



(ln(l + *i))2 



x 1 



S k (xi 



It is more convenient to rearrange this to put the square- 
bracket term on the LHS, and define the quantity Xs ("ex- 
cess speed") by 



Xs(zi,z 2 ) = 



21 ri(zi) 

y(z2) 



(ln(l + 2 1 )) 2 
Z?(l + Zl) 



1/3 



Xl 



2/3 



;(X5) 



where Xs is a ratio of observables, and xi = Dr(z\)/Rc 
as before. (Note one may cancel some powers of Z\ on the 
LHS, but leaving them as above makes both terms in Xs 
well-behaved as z\ — > .) 

This inequality forms the main result of our paper, our 
cosmic speed-trap, which must be obeyed for any chosen 
values Z\ and 2:2 with z\ ^ z 2 , given the following conditions: 

(i) The universe is nearly homogeneous and isotropic with 
a Robertson- Walker metric. 

(ii) The redshift is due to cosmological expansion and c 
is constant. 

(iii) r s is the same comoving length at z\ and 22, and 

(iv) The expansion has never accelerated in the interval 
< z < z 2 . 

If the speed-trap is observationally violated, Xs > 
(xi/S h (xi)) 2/3 at high significance, one or more of assump- 
tions (i)-(iv) above must be false, independent of gravity 
theory or Friedmann equations. To apply this test, we also 
require an upper bound on the RHS, i.e. an upper bound 
on xi for closed models, which we derive below (this is not 
strictly a fifth "assumption", since it follows from observa- 
tional data assuming (i), (ii) and (iv) above). 

In inequality (|15|) . the LHS Xs is formed from a ratio of 
two dimensionless BAO observables d(z\) and 1/(22), while 
the RHS is close to 1 with a weak dependence on curvature: 
the effect of curvature on the low-redshift Dv(zi) is folded 
into the factor containing Sk on the RHS. As noted above, 
this is exactly 1 for flat models and is always < 1 for open 
models, so open models always tighten the speed trap. For 
closed (positively curved) models the Sk factor is > 1, which 
weakens the trap slightly; however, at low redshift zi this is a 
small effect as follows: from the discussion in ijl3.ll for closed 



models we found a conservative lower limit Rc zZ 0.6 c/Hq; 
this leads to X\ ^ ln(l + zi)/0.6, thus for example the RHS 
is ^ 1.013 for 21 = 0.2. The top solid curve in Figure[5]shows 
the resulting upper limit on the RHS of (|15[) assuming the 
very conservative limit Rc 0.6 c/Ho, while the next-to-top 
solid curve shows the limit assuming Rc 1.0c/ffo- 

Thus, if actual observations reveal that Xs > 1.02 with 
good significance, the cosmic speed-trap "flashes" : if so, we 
can then rule out all homogeneous non-accelerated models 
regardless of the detailed expansion history or gravity model. 

In the above Eq. 1151 the square bracket term in Xs is 
given to first order by (1 + §£i)~ . Higher order terms are 
small, and a quadratic approximation is not an improve- 
ment; a slightly better approximation is (1+0.65 21) _1 which 
is accurate to 0.2% for 21 < 0.3. Note that the RHS of OH 
has no dependence on 22; the curvature radius Rc has no 
effect on the observable 2/(22) since y is purely a line-of- 
sight measurement. Therefore, we may choose to measure 
2/(22) anywhere, but if the real universe is accelerating, the 
observed Xs will be maximal when 22 is close to the past 
minimum of 1/(22), at 22 ~ 2 acc . 



3.3 Predictions for ACDM 

In Figure [2] we show predictions for Xs(2i,22) as a func- 
tion of z\ for three ACDM models (dashed) and one wCDM 
model (dash-dot), from substituting Eq. ll3l into ll5l and eval- 
uating H(z) and Dv{z) for the models. For each of these 
plotted curves, 22 is set to 2 acc for that model. The non- 
accelerating upper limit for Xs (the RHS of Eq. ll5p is shown 
as solid lines for several assumed values of curvature radius 
Rc- 

We see from Figure f2] that if the real universe has fol- 
lowed an expansion history H(z) similar to ACDM predic- 
tion, inequality [15] will be violated if 21 is reasonably small 
and 22 is near z a cc ~ 0.75. Essentially, the accelerated ex- 
pansion between 22 and 21 causes the value of H(z)/(1 + 2) 
to be larger at ^ «i than in the past at 22, as in Figure 1; 
this makes Dv(zi) smaller and d(z) larger, compared to any 
non-accelerating model with the same H(z 2 ), so Xs violates 
the limit in Eq. [15] 

As noted above, to maximise the violation we should 
choose 22 to minimise the observed value of y(z 2 ), i.e. the 
redshift 2 acc where H(z)/(1 + 2) had its past minimum; for 
a ACDM model with Q m ~ 0.27, the actual minimum is at 
2 a cc ~ 0.75, but the theoretical y(z) is within 2% of its mini- 
mum value over a rather broad window 0.5 < 2 ^ 1.1: so for 
an observational application of the test, 22 may be what- 
ever is most convenient observationally within this range, 
with only marginal weakening of the trap. 

Turning to the variation of Xs with z\, the predicted 
value of Xs is maximal at 21 = (with a value of 1.185 for 
our reference model C), and slowly declines with z\: thus 
lower 21 is better both to maximise lever-arm in our speed- 
trap, and to minimise curvature uncertainty. However, for 
practical observations z\ cannot be too small since we need 
sufficient cosmic volume to get a robust detection of the 
acoustic feature in the galaxy correlation function £(r) or 
power spectrum P(k); therefore there is a tradeoff between 
Xs which declines with 21, the curvature uncertainty also 
favours smaller 21 , but the available cosmic volume for mea- 
suring d(zi) grows with z\. Thus for an observational appli- 
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Figure 2. This figure shows both sides of inequality I15I as a function of redshift z\. The solid lines show the right-hand side of I15I i.e. 
the upper limit on Xg for non-accelerated models, assuming respective curvature radii Rq = —0.6, —1.0, oo, +1.0, +0.6 in units of c/ Ho 
(bottom to top). The dashed lines show the predicted values of A'.<y(zi , 22) for varying z\ at fixed 22 = z acc , for the same four models as 
in Fig.[T] The three dashed lines show flat ACDM models with Q m = 0.24 (upper), 0.27 (thick), 0.31 (lower). The dot-dashed line shows 
wCDM with U m = 0.32, 10 = -0.85. 



cation of the speed-trap, there is an optimal window around 
0.15 < 0i < 0.35. 

Taking example values Z\ — 0.1, 0.2, 0.3, the concor- 
dance model predicts X s (zi, 0.75) = 1.145, 1.113, 1.088 re- 
spectively. We also note that the value of Xs is fairly sensi- 
tive to the value of fi m : taking example cases from Table [T] 
with Q m — 0.24, 0.27, 0.31 to bracket the plausible range, 
we find that Xs(0,z aC c) = 1.225, 1.185, 1.143 respectively; 
while X s (0.2, z acc ) is 1.142, 1.113, 1.081. For each model, 
Xs — 1 approximately halves from z\ = to z\ ~ 0.27. This 
is because the rate of acceleration grows with time after z aC c, 
so Xs has stronger than linear dependence on qo- 

We note here that the prediction for Xs is independent 
of Ho if all of Q m , Qde, fifc and w are held fixed. How- 
ever, since our example models are approximately CMB- 
matched, a correlation appears, because raising Q m and/or 
w compared to the concordance model requires lowering Ho 
to remain consistent with the CMB; while raising fl m or w 
also leads to weaker acceleration and thus lowers Xs- Thus, 
Xs at a fixed redshift is positively correlated with Ho in 
CMB-matched Friedmann models. 

We also note that for accelerating models Xs remains a 
few percent greater than 1 for the case z\ = 22; this occurs 
because 1/(22) measures the instantaneous expansion rate at 



Z2, while d(zi) depends on the average expansion rate at 
redshifts below zi, which is larger. In principle we could use 
this to test acceleration by measuring d and y from a single 
survey at 21 — Z2, but in practice the curvature uncertainty 
probably disfavours this (see Sec 14.51 for more discussion). 



4 DISCUSSION 

In this section we discuss various aspects of the test above, 
including possible shifts in length r s , useful approximations 
for Dv(z), observational issues, the relation to the Alcock- 
Paczynski ratio and the effect of giant-void models. 

4.1 Possible shifts in r s 

In applying the speed-trap, clearly assumptions (i) and (ii) 
above are very basic; if future observations show the speed- 
trap is observationally violated, we need to be confident that 
assumption (iii) on constancy of r 3 is valid to around ~ 
2%, in order to reject general homogeneous non-accelerating 
models with high confidence. 

We now consider some details which may actually give 
rise to a significant shift in comoving r s between redshifts 
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z\ and Z2\ the main such effects are galaxy bi as, non-linear 
growth of struc ture, redshift-space distortions (|Kaiserll 19871 ; 
iHamiltonl 1 19921 ). and possible effects due to the hump(s) 
sitting on a sloping background power spectrum etc. 

We first note that there is non-negligible evolution in rg 
at high reds hift between Zd ~ 1020 and z ~ 10, as shown by 
Figure 1 of lEisenstein. Seo fc White! 120071 ; the initial BAO 
bump is only in the baryons and photons, and the peak shifts 
slightly as the dark-matter and baryon perturbations align 
together at later times; this implies the late-time BAO peak 
is not exactly at the sound horizon length rs{zd). However, 
after z < 10 the density perturbations in baryons and dark 
matter are very similar. In most real BAO analyses, a matter 
power spectrum from CMBFast or similar is used, together 
with a model for non-linear evolution and an arbitrary linear 
"stretch factor" a, to fit observations; finally, the measure- 
ment is quoted as d(z) = ars(z<i)/Dv{z) where rs and Dv 
are both computed from the reference theoretical model. 
This implies that small errors in the reference model should 
(on average) be absorbed into an opposite shift in a, so the 
final estimate of d(z) should be unbiased. Any shift in the 
BAO length from z ~ Zd to z ~ 10 is included in the refer- 
ence model; therefore, rs(zd) forms essentially a convenient 
fiducial length for intercomparison between models, which is 
close to but not exactly the position of the low-redshift BAO 
peak in the correlation function. For the present work, we are 
only interested in shifts of the BAO scale at z < Z2 ~ 0.75, 
so the above effect cancels. 

Galaxy bias, at least in standard versions, has little ef- 
fect since the BAO scale is very much larger than any scale of 
relevance for galaxy formation; thus bias may affect the over- 
all amplitude of galaxy clustering but cannot significantly 
shift the scale r s . Likewise, non-linear growth of structure 
primarily moves galaxies around on ~ 5 /i -1 Mpc scales; this 
significantly blurs the bump in £(r), and/or erases the higher 
harmonics in the power spectrum, but this is almost sym- 
metrical between inward and outward shifts: the systematic 
shift in the BAO lengthscale is much smaller. 

For the standard model, thes e effects have been inves- 
tigate d fro m both theory by e.g. Eisenstein. Seo fc White! 



|2007l ) and lShoii. Jeong fc KomatsuT 20091 )7 and from large 
N-body simulations bv lSeo et all (120081 ) and lSeo et all (|2010T> ; 



these papers agree that systematic shifts are small, typically 
below the 0.6% level at z = 0.3 and less at higher red- 
shift. They also find that rec onstruction methods based on 
velocity- field reconstruction (|Eisenstein et all 120071 ') can re- 
duce the shift to ~ 0.1%. This will become important for the 
next generati on of ambitious plann ed surveys such as ESA's 
Euclid (e.g. ISamushia et"all Wll ) or NASA's WFIRST, 
which aim to achieve sub-percent precision on BAO observ- 
ables in many redshift bins, but are almost negligible with 
respect to the speed-trap test in this paper. 

We caution that there is a slight level of circular argu- 
ment in the above, in that we are assuming standard cosmol- 
ogy to limit the shift in r s , and then using this to reject non- 
standard non-accelerating models; it remains possible that 
a model with non-standard gravity could produce a much 
larger shift in r 3 than the standard cosmology. However, 
non-standard models producing a gross ~ 10% shift in r s 
since Z2 ~ 0.75 would almost certainly produce large levels 
of redshift-space distortion, and give strong inconsistencies 
between the angular and radial measurements of r 3 at low 



redshift. If both the redshift-space distortion pattern and 
the radial and angular measurements of r s are measured to 
be consistent with standard ACDM, this would strongly sug- 
gest that the true shifts in r s should not be much larger than 
the percent level effects predicted by the standard model. 

4.2 Approximations for Dv 

As an aside, we also note that in nearly-flat CDM-like 
models, an accurate approximation to Dv(z) at moderate 
redshift is given by Taylor-expanding 1/H(z) around z/2 
(rather than zero), and substituting in the integral Eq. 2] 
this makes z 2 terms vanish, and leads to the approximation 



D v (z) 



[H 2 (z/2)H(z)] 1 / 3 



+ 0(z s 



(16) 



in practice the first term is surprisingly accurate for ACDM 
models, with errors < 0.1% compared to the numerical result 
for z < 0.5. (See Appendix A for evaluation of the third- 
order term, and explanation why it is small). 
A simpler approximation is 



D v (z) 



H(2z/3) 



(17) 



this is slightly less accurate than the previous approxima- 
tion, but still accurate to < 0.4% for z < 0.5, better than the 
mid-term precision on observables. (For open zero-A models 
these approximations are less good, with errors up to 2%). 

While it is straightforward to evaluate Dv and Xs nu- 
merically for any given model, the main value of this approx- 
imation is that it tells us that a measurement of z d(z) at 
low redshift is quite close to a measurement of r 3 H(2z/3)/c 
; substituting this into (115|) . along with the approximation 
(1 + fzi)" 1 for the square-bracket term, gives simply 



Xs(zi,z 2 ) ; 



+ (5*1 J 



Zl 



H(z 2 ) 1 + 



a(z 2 ) 



(18) 



and inequality [3] tells us this should be less than 1 for non- 
accelerated models. Unlike our upper limit Eq. (|15[) this ex- 
pression is not rigorous, but this gives a simple and fairly 
accurate approximation for what Xs is measuring, i.e. it 
is closely related to the ratio of expansion rates a at |«i 
compared to z 2 . 



4.3 Observational advantages 

One possible objection to this test is that it is comparing 
two related but slightly different observables, i.e. a spherical- 
average scale at z\ with a radial scale at z 2 . Why have we 
done this, rather than comparing two measures of d{z) or 
two measures of y(z) at two different redshifts ? 

It is well known that comparing y(z) at two differ- 
ent redshifts provides another direct test of acceleration. 
The main difficulty is observational, since for our baseline 
model, y(zi) only grows to l.ly(z a cc) at rather low redshift 
Zi ~ 0.16. Furthermore, for a given survey, a radial-only 
measurement of r s has a statistical error roughly y3 worse 
than a spherical average measure. Even if we had a 3n stera- 
dian redshift survey at Z\ ~ 0.16, we may not do much better 
than 3% statistical error on y, a 3a violation, and we would 
like to get above 5a for a decisive result. Using d(zi) instead 
of y gives two substantial advantages: firstly d(z\) effectively 
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measures H at ~ 2zi/3, giving more lever-arm on the low- 
redshift acceleration; so a measurement of d(z± = 0.24) is 
similar in content to a measurement of y(zi = 0.16). Sec- 
ondly there is the obvious gain that d uses 3 spatial direc- 
tions instead of 1. Thus for a fixed thickness of survey shell, 
the former measure has around 9/4 times more available vol- 
ume and 3 independent axes, s o the cosmic variance limit 
should improve by a factor ~ -y/27/4 ~ 2.6, which is a very 
important practical advantage. 

In contrast, comparing d(z) at two different redshifts 
suffers from potential major uncertainty in cosmic curva- 
ture at the high redshift 22. At 22 ~ 0.75 there is ample 
available volume for a precision meas urement of y, and am - 
bitious future probes such as Euclid (|Samushia et all l201lh 
plan to push to statistical errors < 0.75% on both y and d, 
in each of many bins of width 0.1 in redshift. Thus, at Z2 
the cosmic variance is minimal for a wide-area survey, so the 
radial measure is preferable because it is independent of the 
curvature nuisance parameter. Also, d(z2) depends on the 
full history of H(z) back to Z2, which complicates the issue 
of deriving an inequality. 

In our proposed comparison, we have constructed a ra- 
tio Xs using d(zi) at low redshift and y{z2) at the higher 
redshift, to circumvent both of these problems: the potential 
cosmic-variance limits are probably around 1% on d(0.24) 
and significantly less on 2/(0.75), so this test can (given am- 
ple data) deliver a standalone rejection of homogeneous non- 
accelerating models at ~ 7<r significance level. This can be 
further improved by using several independent redshift bins, 
e.g. 2i = 0.15, 0.25 and z 2 = 0.65, 0.75. 



4.4 Future Observations 

As noted above, there already exist measurements of th e 
numerator on the left of Eg. 1141 from IPercival et all (2010); 
they quote values of d(0.2) = 0.1905 and d(0.35) = 0.1097, 
with approximately 3.3% error on each. For the numerator 
z d(z) in inequality |15]) these give 0.2 d(0.2) = 0.0381 and 
0.35d(0.35) = 0.0384. 

As yet there is no available measurement of radial BAOs 
at 2 > 0.5 with which to actually calibrate our speed-trap, 
but these are expecte d soon \^ from the recently completed 
AAT WiggleZ survey (|Blake et alllgOlOh. and in a few years 
from the ongoing BOSS survey ( White et all l201lT ). It is 
currently unclear whether the final WiggleZ survey covers 
enough volume to separately measure the radial component 
as required here, but BOSS should very likely achieve this; 
the upper redshift limit of BOSS is ~ 0.65, so this is close 
enough to 2 acc to be useful. 

For ACDM, the predicted value of 4/(22) near its min- 
imum is approximately 0.0302 for Q m = 0.27 and Ho = 
70kms _1 Mpc -1 . For reasonable variations of parameters, 
we now show that if we assume a flat universe then 2/(22) is 



well constrained by CMB observations: it is well known that 
for flat models with varying Q. m , h, w there is a tight corre- 
lation between th e age of the Univ erse, to, and the CMB 
acoustic scale £a (|Knox et alH200lh ; and it turns out that 
there is also a tight correlation between these and the value 
of H at intermediate redshift, with a pivot point occurring at 
2 ~ 0.8 (see Figure [1]). This is partly a coincidence, because 
for moderate parameter variations around the concordance 
model, to scales tx Sl~°' 3 /i _1 , while £a scales as Q.^' 15 h~ ' 5 . 
For the value of H(z), we note that as 2 — > this scales as 
h independent of Q, m , w, while at 2 > 3 where dark energy 
is negligible, H(z) scales oc 0,^'h 1 . Therefore, there exists 
a pivot-point at intermediate redshift where H(z) scales as 
Q^h (i.e. inversely to to), and this pivot redshift turns out 
to be 2 ~ 0.85 for A models. For w > — 1 the pivot red- 
shift is somewhat lower, but for near-flat Friedmann mod- 
els the value H{z — 0.75) is better constrained by WMAP 
data than the local Ho; and fixing to ~ 13.75 Gyr constrains 
H(z = 0.75)/(1.75) = 59.2 km s" 1 Mpc" 1 within « 0.8%, 
which in turn leads to a tight prediction for y(0.75). 

(As an aside, there is a corollary that if some future 
method could give a direct measurement of H{z = 0.75) in- 
dependently of rs, this would produce another strong con- 
sistency test of standard ACDM. This may be possible in 
principle using methods such as differential-age measure- 
ments of early-type galaxies, or lensing measurements with 
source and lens close in redshift, but this will require a major 
advance in precision over current data). 

Assuming some future y measurement turns out at the 
concordance value 2/(22) ~ 0.0302, we would then obtain 
measurements Xs ~ 1.12 and 1.05 at z\ — 0.2, 0.35 respec- 
tively. The error on y(z2) must be added in quadrature to 
the current 3.3% error on d(zi), but if the former is around 
2% then we can anticipate a fairly clear violation from the 
2 = 0.2 value, and a somewhat less significant violation at 
2 = 0.35. 

The prospects are good for improving on the current 
resul ts: the projections for the BOSS survey (Whi te~et all 
l201l|) are for 1% precision on d(z = 0.35), and precision 
of 1.7% on 7/(0.6). Adding the above errors in quadrature 
leads to around 2% precision on Xs, with a predicted value 
~ 1.077, thus nearly a 4<r proof of acceleration. BOSS may 
also do better using the larger value of Xs at z\ ~ 0.2, but 
projected precision on d(0.2) is not quoted separately. 

Next-generation surveys in the planning stage such 
as BigBOSS, Euclid or WFIRST should substantially im- 
prove on the higher-redshift measurement, reaching sub- 
percent precision on y(z2)- The low-redshift d(z\) measure- 
ment is ultimately limited by cosmic variance, but extend- 
ing the BOSS survey to the Southern hemisphere can give a 
straightforward improvement by a factor of \pi, or probably 
more if denser sampling of galaxies is used. Further improve- 
ments are possible in principle using HI or near-infrared se- 
lected surveys which can cover > 80% of the whole sky, 
compared to ~ 50% for visible-selected surveys. 



1 Soon after the submission of the first version of this paper, 
three new measureme nts of the BAO feat ure appeared: one from 
6dFGS at z ~ 0.1 inlBeutler et~al i201lT ), one from WiggleZ at 
z = 0.6 inlBlake et aJ 1 2011 ). and one from SDSS photo-z's at 2 = 
0.55 in lCarnero et all 201S ). All of these show good consistency 
with the concordance model, but do not yet measure the radial 
component as required here. 



4.5 Comparison with the Alcock-Paczynski test 

We note here that our ratio Xs may be considered as a 
genera lised version of the classic test of lAlcock fc Paczvnskil 
( 1979), hereafter AP: the AP ratio was defined to be Rap = 
Az/zA6, which in our notation becomes 
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Rap(z) 



(1 + z)D a {z)H(z) 



(19) 



If we choose zi = Z2 in Eq. [13] above and substitute Eq. [T] 
for Dv, we then obtain 



(20) 



Zld(zi) , . , ,_ 2 /3 

— -. = (1 + zi)Rap(zi) 

y{zi) 

thus Xs(z\,z\) contains the same information as Rap(zi) 
combined with a function of Z\\ substituting the above into 
Eq. [15] gives a lower limit on Rap for non-accelerating mod- 
els, which is 



Rap(zi) > 



(l + zi)ln(l + Zi) S k ( Xl ) 



(21) 



Zi X\ 

It is well known that if we assume the Friedmann equa- 
tions, the AP test at high redshift provides a strong test 
for A or dark energy: however, if we drop the Friedmann 
connection between curvature and matter content, then at 
z > 0.5 the AP test becomes mostly degenerate between 
acceleration and curvature. At lower z < 0.4, we may use 
the approximation (1 + z)Da{z) ~ cz/H(z/2) from above, 
which leads to Rap(zi) « H(zi)/H(^-). This does have 
more sensitivity to acceleration than curvature, but is not 
ideal for the following reason: at small z\ the AP ratio suf- 
fers from a short redshift lever-arm, while at z\ > 0.4 the 
ratio mainly probes the regime of sluggish acceleration at 
z > 0.2. The AP ratio at Z\ ~ 0.4 may provide a useful test, 
but will probably require sub-percent level precision on both 
observables to get a decisive result. 

Compared to the AP test, the use of two widely-spaced 
redshifts in Xs requires the added assumption that rs has 
minimal evolution between Z2 and z\, but enables a much 
longer effective time lever-arm, giving a larger acceleration 
signal while keeping the curvature sensitivity very small. 

4.6 Inhomogeneous Void Models 

Recently there has been some interest in models which pro- 
duce apparent acceleration without dark energy, by plac- 
ing us near the centre of a giant underdense spherical 
vo id, with a Lem aitre-Tolman-Bondi metric; examples are 
in lTomital (2009) and references therein. These models have 
several problems such as severe fine-tuning of our loca- 
tion very close to the void centre, and probable inconsis- 
tency with limits on th e kinetic Sunyaev-Zeldovich effect 
( Zhang & Stcb binsl 1201 lh ; however it is interesting to note 
how Xs behaves in such models. A recent confrontation of 
giant-void models with BA O observables has been done by 
iMoss. Zibin fc Scottl (|201lh : they find that void models with 
profiles adjusted to match SNe and CMB observations have 
a Az« which is > 30% smaller at z ~ 0.5 — 0.7 compared to 
ACDM. Those specific cases would have X s (0.2, 0.75) > 1.4, 
which is substantially larg er than any reason able dark- 
energy model; thus, IMoss. Zibin fc Scottl (|201ll ) show that 
giant-void models matched to angular distances and the 
CMB appear to suffer from severe "overkill" in radial BAO 
measurements. 

The parameter space of possible void models is very 
large, so other void models may look more similar to ACDM , 
but we note that the test of lClarkson. Bassett fc Lul {2008) 
can be used to test for homogeneity without assuming GR. 
They show that if we have both angular and radial BAO 



measurements spanning a range of redshift, there is a con- 
sistency relation which must be satisfied by homogeneous 
models but is usually violated by giant-void models. Thus, 
assumption (i) above becomes observationally testable us- 
ing future BAO observations, though this probably requires 
observations spanning more redshifts than the Xs test here. 



5 CONCLUSIONS 

We have proposed a new and simple smoking-gun test for 
cosmic acceleration using only a comparison of the baryon 
acoustic oscillation feature at two distinct redshifts ~ 0.2 
and ~ 0.75. The main result of our paper is inequality ()15|) 
relating the two dimensionless BAO observables, which must 
be satisfied for any homogeneous non-accelerating model, 
but will be observationally violated by w 10% in models 
with an expansion history close to standard ACDM. 

Clearly, our proposed measurement has advantages and 
disadvantages: the main advantages are extreme simplic- 
ity and model-independence, i.e. if the inequality (|15|) is 
violated, we can rule out essentially all homogeneous non- 
accelerating models in one shot, without assuming any par- 
ticular gravity theory or parametric form of H(z), and in- 
dependent of supernova and CMB observations. 

The main drawback of our test is that it is essentially 
one-sided: if inequality [TS] is observationally violated, we 
have proved (given some basic assumptions) that acceler- 
ation has occurred during < z < Z2 and have a rough 
quantification of the amount, but no more details about the 
underlying cause or the details of the expansion history. 

If we assume GR and the Friedmann equations hold, 
and that r s has the value which is accurately predicted from 
CMB analysis, then we have much more statistical power: 
future measurements of BAOs in many redshift bins may 
be used to reconstruct the detailed form of the expansion 
history H(z); this can be integrated to give predictions of 
Dr(z), and comparison with the measured transverse BAO 
scale giving Da(z) can constrain spatial curvature inde- 
pendent of the CMB; while comparison of Da with Dl(z) 
from SNe can check the distance-duality or Tolman relation 
(Dl/Da) 2 = (1 + z) 4 . All of this can give much more power- 
ful cross-checks and parameter estimates than our simplified 
one-sided test. 

However, our proposed cosmic speed-trap seems to pro- 
vide a valuable addition to the set of cosmological measure- 
ments, due to its bare minimum of assumptions. This pro- 
vides a strong motivation for future improved BAO mea- 
surements specifically near redshifts ~ 0.25 and 0.75; this 
should preferably include a low-redshift survey comparable 
or superior to BOSS in the Southern hemisphere to minimise 
the cosmic variance in the local measurement. 
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APPENDIX A: THE APPROXIMATION FOR 

D v 

We here add a note which explains the surprisingly good 
accuracy of approximation[l6]for Dv (z) at fairly low redshift 
z < 0.4. As noted, in the integral Eq. [4]for Dr, it is helpful 
to Taylor-expand the function 1/H(z) around the mid-point 
of the integral at zi/2, then integrate: this naturally makes 
terms with odd- integer derivatives otl/H integrate to zero, 
and leads to 



D R {zi) 



(Al) 



where prime denotes d/dz. We now need the second deriva- 
tive (1/H)" evaluated at 4^. Defining the usual deceleration 
parameter q and the jerk parameter j fe.g. lAlam et alll2003l ) 

as 



d 2 a/dt 2 
aH 2 



d 3 a/dt 3 
aH 3 



we can rearrange these in terms of d/dz to get 



d 2 H 

dz 2 



dH H 
-dl = — {1 + q) 
Using these we obtain 

d 2 fj_^ _~j + 2 + 4q + 3q 2 

dz 2 



H 



r0' 



(A2) 



(A3) 



(A4) 



(l + z) 2 H 

For the case of flat ACDM models, j = +1 indepen- 
dent of parameters (assuming radiation density is negligible) 
|Rapetti et all [2007). thus the numerator in Eg . I A4I has ze- 
ros at q = — 1/3 and q = — 1. For f2 m near the concordance 
model, q passed through —1/3 in the fairly recent past at 
2 ~ 0.3, so the numerator is significantly smaller than unity 
at low redshift. This explains qualitatively the very good 
accuracy of approximation \W\ near the concordance model, 
even up to significant redshifts z ~ 0.5. 

This paper has been typeset from a TgX/ DTgX file prepared 
by the author. 



© 2011 RAS, MNRAS OOO.riHTol 



